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Abstract 



Aims. On the basis of the PPMXL star catalogue we performed a survey of star clusters in the second quadrant of the Milky Way. 
Methods. From the PPMXL catalogue of positions and proper motions we took the subset of stars with near-infrared photometry 
from 2MASS and added the remaining 2MASS stars without proper motions (called 2MAst, i.e. 2MASS with astrometry). We 
developed a data-processing pipeline including interactive human control of a standardised set of multi-dimensional diagrams to 

■ determine kinematic and photometric membership probabilities for stars in a cluster region. The pipeline simultaneously produced 
_ the astrophysical parameters of a cluster. From literature we compiled a target list of presently known open and globular clusters, 

\^ ■ cluster candidates, associations, and moving groups. From established member stars we derived spatial parameters (coordinates of 

centres and radii of the main morphological parts of clusters) and cluster kinematics (average proper motions and sometimes radial 
velocities). For distance, reddening, and age determination we used specific sets of theoretical isochrones. Tidal parameters were 
obtained by a fit of three-parameter King profiles to the observed density distributions of members. 
^ . Results. We investigated all 871 objects in the 2nd Galactic quadrant, of which we successfully treated 642 open clusters, 2 globular 

5_( ■ clusters, and 8 stellar associations. The remaining 219 objects (24%) were recognised by us to be nonexistent clusters, duplicate 

■ entries, or clusters too faint for 2MAst. We found that our sample is complete in the 2nd quadrant up to a distance of 2 kpc, where the 
\ average surface density is 94 clusters per kpc^. Compared with literature values we found good agreement in spatial and kinematic 

data, as well as for optical distances and reddening. Small, but systematic offsets were detected in the age determination. 

Key words. Galaxy: globular clusters: general - Galaxy: open clusters and associations: general - Galaxy: stellar content - Galaxies: 
fundamental parameters - Galaxies: photometry - Galaxies: star clusters 
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[ 1 . Introduction cannot be resolved), and we have to rely on more sophisticated 



models such as simple stellar population (SSP) models to esti- 



1^ . Star clusters are important representatives of the galactic popula- j^ate fundamental parameters of the clusters. 
j3 ■ tion. Their fundamental astrophysical parameters, such as ages j„ ^^e Galaxy, stai- clusters are resolved into individual 
(N| . and masses, can be determmed much more reliably and accu- ^^^^^^ the bulk asti-ophysical parameters can be derived 
^ . rately than for individual stars. In the past, the efforts of as- properly. :„ the case of open clusters we have to pay for 
•' ; tronomers concentrated on the analyses of individual clusters j^is huge advantage with a limited horizon of a few kilopar- 
IIh . and their properties in our Galaxy and the Magellanic Clouds. ^^^^ ^q^^^ the Sun. In the past, much effort has been de- 
^ ■ During the last few decades the interest has shifted to studies ^^^^^ ^^Ueet different (usually heterogeneous) data on in- 
^ : of large samples of clusters. Galactic and extragalactic, which ^^^^^^^1 clusters to p r ovide a basis for study i ng their pop- 
. P_ . yielded clues for our understanding of the history of star for- ^j^jj^n (B^ [IM iBecker & Fenkartl [iM [S 
mation and the evolution of galactic populations over large spa- ^^^^^ ^^gg) ^he most recent collections of this nat ure are 
tial scales. Unlike in our own Galaxy, the distribution of open ^^e updated on-line list of cluster data of iDias etalJ (l2002l) 
star clusters in nearby galaxies can be observed over the whole j^e database on Galactic clusters WEBDAQ. Nowadays, 
body of a galaxy, which permits direct access to the spatial dis- j^rge all-sky photometric and/or kinematic surveys have opened 
tnbution within the host galaxy. As examples we mention re- ^^e possibility to derive homogeneous data sets of the pa- 
cent studies of various properties of extragalactic clusters in dif- rameters of Ga lactic clusters. Especially the 2MASS survey 
ferent environments that involved hundreds to a few thousands jskrutskie et al.1 120061) has considerably contributed to the im- 
of objects (see e.g. iCaldwelLet al. 2009; Cantiello et al.ll200| provement of our knowledge on different aspects of Galactic 
[Moraetal. 2009; Vansevicius et al. 2009; Harris et al. 2010t clusters. Among o t her studies, w e men tion the w ork by 
San Roman et al. . ,201^ , Werchan & Zaritsky . .2011). The price 'Dutra & Bica ( 2000), Dutra & Bica (2001), Dutra et al. (^Oo|^ 



for this advantage is a poorer resolution (mostly individual stars frpebrich et al. (2007), Froebrichet al., (.201Q) . .Glushkova et : 
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(1201 Oh . iBukowiecki et alJ (l201lh . and iTadrossI (l201lh . which 
deal with large samples of the Galactic star clusters whose deter- 
mination of the cluster parameters is based on 2M ASS data. We 
also mention studies of smaller samples of clus ters dlvanov et al.l 
I2OOIIK0POSOV et al.ll200ilMercer et alj|2005l) . 

Obviously, the reliability of the results from photometric 
data can be considerably improved if kinematic information is 
also included in the membership determination. However, an in- 
creasing number of clusters with homogeneous parameters alone 
is not sufficient for studies of the cluster population: knowing the 
statistical properties of cluster samples (how complete, how rep- 
resentative) is another basic requirement. 

In our previous study we used the ASCC-2.5 (I Kharchenko! 

I2OOI catalogue to identify known clusters and systematically 
search for new ones. This resulted in a sample of 650 local open 
clusters with individual membership information for cluster stars 
based on kinematic, photometric, and spatial criteria. For each 
cluster a set of homogeneous parameters was determined in- 
cluding spatial (positions, distances), structural (apparent sizes), 
kinematic (proper motions and sometimes radial velocities), 
photometric (reddenings, integrated magnitudes and colours), 
dynamical (King' ti dal parameters and masses ) a nd evolutionary 
(ages) parameters ( Kharchenko et al. 2005a||g; Piskunov et alj 



120071; iPiskunov et al.ll2008l: lKharchenko et al. 2009). The results 
are summarised in the Catalogue of Open Cluster Data (COCD), 
available at the CDS online archive. The cluster sample is com- 
plete within about 0.85 kpc and describes the local population 
of open clusters. However, owing to the bright limiting magni- 
tude of ASCC-2.5, this is sufficient only for studying the local 
inter-arm cluster population of the Galactic disk. To reach larger 
distances one needs a deeper survey or a shift to less obscured 
wavelengths. 

The current project "Milky Way Star Clusters", hereafter 
referred t o as MWSC, is b ased on the new all-sky catalogue 
PPMXL jRoser et al. 112010'), which goes considerably deeper 
than ASCC-2.5, and provides a substantial expansion of the ac- 
cessible volume. We aim at identifying different cluster-like ob- 
jects (open and globular clusters and candidates, associations, 
moving groups, cluster remnants, etc. from a candidate list) in 
the PPMXL and at determining their fundamental parameters. 
We use an automated pipeline for membership and cluster pa- 
rameter determination which is, additionally, under human con- 
trol at each stage. We expect to considerably increase the number 
of clusters having a homogeneous set of cluster parameters and 
reach the nearby spiral arms with the completed part of the sur- 
vey. The data products will include information on membership 
of stars in cluster areas (catalogue of stars, MWSCS), and for 
each confirmed cluster we will provide a set of homogeneous 
basic cluster parameters (MWSCD), together with an atlas of 
cluster diagrams (MWSCA). We started at I = 90° and are mov- 
ing along increasing galactic longitude. 

The current paper releases the data on all clusters identified 
within the 2nd Galactic quadrant (I = 90° . . . 180°). In Sect. |2] 
we specify the basic input data. The pipeline of the cluster pa- 
rameter determination is described in detail in Sect. [3] In Sect.|4] 
we estimate the statistical properties of the cluster sample so far 
and give preliminary results on the spatial ditribution of open 
clusters in this sector of the Milky Way. Sect. |5] summarises the 
first results of the project. 

2. Data 

For our study, t he basic data we re taken from t he all-sky cata- 
logues PPMXL (lRoseretalJl201Q) and 2MASS dSkrutskie et alJ 



I2006h . PPMXL gives positions and proper motions in the 
International Celestial Reference System (ICRS) and low- 
accuracy photometry from USNO-Bl.O (Monet et al. 2003) for 
about 900 million objects down to V ^ 20. For some 400 mil- 
lion entries the catal ogue contains accura te J, H, Kg magni- 
tudes from 2MASS dSkrutskie et all l2006h . We note that the 
2MASS limiting magnitudes vary significantly due to crowding. 
This subset was used to verify the cluster and determine clus- 
ter parameters in the astrometric and photometric systems (of 
PPMXL) that are homogeneous over the whole sky. 

2.1. Construction of 2MAst 

For about 70 million stars from 2MASS no counterparts are 
given in PPMXL. These stars are either severely reddened along 
the Galactic disk and bulge, or are extremeley red and too faint 
in the optical to be detected in USNO-B 1 .0. Another portion of 
stars is missing in PPMXL because they are blanketed by large 
images of bright stars with ghosts and spikes appearing on the 
Schmidt plates. To avoid incompleteness of the survey in such 
areas, we kept these 2MASS stars with their 2MASS photome- 
try and used them for photometric selection of cluster members. 
In case of extended objects, PPMXL and 2MASS (though the 
latter less often) may include spurious entries and artefacts, and 
about 12% of 2MASS objects have more than one USNO-Bl.O 
counterpart ( Roser et al.ll2010l) . This occurs mainly because of 
the USNO-B 1 .0 matching problems in areas where the Schmidt 
plates of the Palomar survey overlap. In these cases we averaged 
the corresponding proper motions of the PPMXL. 

The resulting data set 2MAst (stands for 2MASS with astro- 
metric data) includes about 47 1 million stars. For about 399 mil- 
lion stars it contains the coordinates and proper motions from 
PPMXL, and near infra-red (NIR) photometry and flags from 
2MASS. For the remaining stars, only the original information 
from 2MASS is available. The limiting magnitudes are 17.1, 
16.4 and 15.3 mag in J (1.25 nm), H (1.65 nm), and Ks (2.17 
nm), respectively. Typical rms errors in proper motions range 
from better than 2 mas/yr for the brightest stars with Tycho-2 
(lH0g et alj l2000h data to 10 mas/yr or even 1 5 mas/yr in the 
region soutii of S < -30° (iRoser et al.ll2010h . For the photo- 
metric data, the la uncertainty is typically better than 0.03 mag 
(Skrutskie et al. 2006) for bright stars (< 13 mag) and increases 
up to 0. 1 mag for fainter stars. 

2.2. Clusters and stars 

Although the main targets of our study are Galactic open clus- 
ters, we also considered such representatives of the Galactic pop- 
ulation as stellar associations, embedded clusters, and globular 
clusters, which can also be identified in 2MAst. 



2.2.1 . Compiling tine cluster list 

The target list of 3784 entries with initial cluster parameters was 
compiled from sources available in the literatu re. As the primary 
source we used the data from the COCD dKharch enko et al] 
2005a b), which was given the highest pr iority. For more o pti- 
cal clusters the data were taken from the Dias et all d2002 l) list 
(Vers ion 3.1, 24/nov/2010). Kn own associations were retrieved 
from iMelnik & DambisI d2009l) . For clusters detected in the 
NIR the information came fromlBic^t al. (2003b). Dut ra et al.l 



I2003h. iF roebrich et al 



IBukowiecki et al.i 6201 1 



d2007h . IFroebrich et al] (|20ll), and 
). Globular clusters were selected from 
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Table 1. Statistics on available data for star clusters in the Milky 
Way 



Parameter 


Number 


Number of entries 


3784 


Angular size 


3648 


Proper motion 


957 


Radial velocity 


730 


Distance 


2175 


Reddening 


1724 


Age 


1562 


King tidal parameters 


1830 



the catalogue bv lHarrisI dl 9961) [edition 2010 R Additionally, we 
incorporated data on radial velocities from iKharchenko et alJ 
(2007) as well as s uppl ementary data on em bedded clusters from 
fBica et al] (l2003a) and Lada & Lada (2003) and on stars in asso- 
ciations from Humphreys (1978). The statistics of the available 
data is given in Table [1] We note that the cluster parameters in 
Table[T]are highly heterogeneous because they are based on dif- 
ferent kinds of observations as well as on different methods of 
membership defintion and parameter determination. 

The compiled list of targets was sorted with increasing right 
ascension and numbered from MWSC 1 to MWSC 3784. We 
kept this identification even if the coordinates of cluster centres 
turned out to be changed in the process of the study. 

2.2.2. Selection of areas and stars 

Around each cluster we initially defined a circular area with a 
radius of 

Ta = Tel + Tadd , 

where Vd is taken from the literature, and Tadd = 0?3. In a few 
cases, where Vd is unknown, we adopted = Vadd as a first 
approximation. If rd and/or the coordinates of the cluster centre 
changed during the data processing, we repeated this step with 
new input values. 

In the areas defined above we selected in 2MAst only those 
stars with flags Rflg (the 2nd triple of the flags in 2MASS) set 
to 1, 2, or 3 in each band, i.e. the stars with the best-quality 
detections in photometric and astrometric data. For a few stars 
spectral types are available from ASCC-2.5. These data were 
also incorporated since this additional information can be useful 
for the determination of cluster parameters. 



3. Construction of the pipeline for cluster 
parameter determination 

The main task of the pipeline is separating of cluster stars from 
the field by use of kinematic, photometric, and spatial criteria 
and determining cluster parameters from the data of the most 
probable cluster members. The pipeline is conceived to work it- 
eratively, a nd, in principle, its philosophy is very similar to what 
we used in IKharchenko et akl (12005 a). However, a switch from 
the optical (B, V) photometry to the NIR (J, H, Kg) photometry 
required considerable modifications of the pipeline. 



3. 1. Isochrones and reference sequences in 
colour-magnitude diagrams (CMD) 

The NIR isochrones for the 2MASS photometric system were 
_computed iising the Padova web-server CMD2.fl based on the 
Vet all (2008h calculations for Z = 0.019. To improve 



the agreement of the lower part of the ZAMS with observations, 
we substituted the isochrones in the dom ain of late-K— M-stars 
by recent calculations of the Pisa group dPa Rioll20 10': Moronjl 
[2010), taking also into account the latest Allard modelfl of M- 
dwarf atmospheres. The pre-main sequence (Pre-MS) portion of 
the isochrones for logt < 8.0 was interpola ted from Pre-MS 
tracks provided by Lionel Siess' web-serve fl dSiess et alJl2000h 
for Z = 0.02. The derived \ogg, logL/L© and T^/f were then 
transformed to the JiJ A's-system with help of the Padova trans- 
formation tables. 

We also used theoretical isochrones to define the zero- 
age main sequence (ZAMS) and terminal-age main sequence 
(TAMS) sequences in the CMD. For the brighter part (earlier 
than K-type), the ZAMS is represented by the Padova isochrone 
for t = 0.001 Myr, whereas for the fainter part, the Pisa Pre- 
MS isochrone for t = 650 Myr and {Y, Z) = (0.28, 0.016) was 
chosen. The latter pro vides a better agr eement with observations 
of Hyades M-dwarfs (iRoser et al. ''201 1). Thus, at low luminosi- 
ties, the ZAMS defines a lower envelope of the cluster MS and 
coincides with the locus of evolved Pre-MS stars at the age of 
650 Myr To avoid a discontinuity due to the use of different 
models, we coiTected the Pisa isochrone, assuming that both se- 
quences must demonstrate the same mass-luminosity relation in 
the log m/mQ ~ —0.1 . . . 0.05 mass range. The respective cor- 
rections are quite small, AA/^, = 0.11, A(J - Ks) = 0.018. 
As a result we have a smooth ZAMS. The TAMS sequence is 
constructed from models of stellar evolution dating the end of 
the hydrogen-burning phase in the cores. Both curves are shown 
in the Kg, {J — Kg) diagram in Fig. lC.2l in Appendix C. 

3.2. Interstellar extinction 

The three-band photometry from 2MASS allows one to deter- 
mine the parameters of interstellar extinction from the photo- 
metric data only. The coiTesponding r elati ons were taken from 
the average extinction law of ICardelli et al.l dl989i) . whereas 
the transformation of ii'-magnitudes to Kg was taken from 
lDutra&Bical(l200lh : 

Ak/Ak^ = 0.95, 

Av/E{B~V) = 3.1, 

E{J - Kg)/E{B -V) = 0.480, (1) 

AKjE{J~Kg) = 0.670. 

The reddening-free parameter Qjhk^ was adopted in the form 

QjHK^={J -H)-x{H -Kg), (2) 

where the slope of the reddenning line x is in first approximation 
equal to the ratio of the colour excesses 



X^E{J- H)/E{H - Kg 



(3) 



Apart from parameters of the photometric system, x depends 
on the stellar spectral type, on properties of the absorbing mat- 
ter along the line of sight (i.e. on the extinction law), and on 



http://www.physics.mcmaster.ca/resources/globular.html 



^ http://stev.oapd.inaf.it/cgi-bin/cmd 

^ http://perso.ens-lyon.fr/france.allard/ 

' http://astropcO.ulb.ac.be/ siess/database.html 
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the extinction value itself. Therefore, it is desirable to derive 
the slope directly from observations whenever possible. Based 
on earlier publications, Mathis (199Q) showed that the ratio 
E{J - H)/E{H - K ') vari es betwe en 1.61 and 2 09, but pref- 
ered the former value. IStraizvs & Laugalvsl (|2008|) have found 
from observations of red-clump giants in the inner Galaxy that 
for extinction Ay < 12 the slope is close to x—^- Before ac- 
cepting this value, we tested the effect of variations of x and 
found that only x = 2 provides independence of Qjhk^ from 
the reddening. Lower values lead to an increase of Qjhk^ with 
E{J — Kg), while higher values decrease Qjhk^ ■ 

Before applying the above relations we tested them with help 
of available photoelectic UBV-photometry in open clusters lo- 
cated in different directions of the Milky Way and found no bias 
between our determinations and observations. The subsequent 
comparison with literature data carried out in Sect l3.4.3l has con- 
firmed the correctness of this procedure. 

3.3. Membership probabilities 

Following the app roach applied earlier for the COCD clusters 
dKharchenko et al. 2004), we computed both kinematic and pho- 
tometric membership probabilities for stars within the cluster ar- 
eas. They were determined by taking into account the accuracy 
of the proper motions and of the stellar magnitudes. 

We define the kinematic probability P^^^ of the «th star to 
belong to a cluster as 



exp ■ 



My My 



, (4) 



where /i^ y and 'p.^ y are the components of the proper motion of 
the ith star and the mean proper motion of the cluster, respec- 
tively. The parameter e^i was set to 1.5 mas/yr for each star that 
had mean errors of proper motion smaller than this threshold in 
2MAst. For all other stars, e^i corresponds to the mean proper 
motion errors given in 2MAst. 

The photometric probabilities Pjh and Pjk are computed 
in two CMDs, Kg, [J - H) and Kg, (J - Kg), respectively. 
Aga in, we followed the approach developed in lKharchenko et all 
(1200 4) for the photometric selection of cluster members. For 
each isochrone we defined the most probable location of clus- 
ter members as an area in between the corresponding smoothed 
isochrones for single stars and for unresolved binaries of equal 
mass. It follows that at a given magnitude Kg there exist two val- 
ues of colour Cfe (blue) and c.^ (red), which limit this domain in a 
CMD. A star i with a colour d is assumed to be a photometric 
member with a probability P^ — 100% if Cb < < Cr- For the 
remaining stars we computed the photometric probability as 



P' 



cxp • 





'Ac'' 








1 



(5) 



where Ac* — & — Cb for < Cb, and Ac' = c* — c^ for c* > c^. 
The quantity e* is the mean error of the colour index c*. It is 
calculated from the individual rms errors of the magnitudes. A 
minimum rms error of a magnitude of 0.03 mag is assumed. In 
this way we obtained two photometric probabilities, Pjh and 
Pjk for each star 

Finally, we defined a combined probability P that takes into 
account all aspects of the membership selection procedure, i.e. 
kinematic, photometric, and spatial selection criteria {Pg) 



where the factor P, is set to 1 within the cluster radius 
r2, and to otherwise (see Sect. 3.4.1). The probabilties 
Pkin , Pjh , Pjk , Ps are given in the catalogue of cluster stars 
MWSCS, and the user may feel free to define his/her own selec- 
tion criteria^ 

As in iKharchenko et al.l (l2004l) . we classified the quality 
of cluster membership. We considered three classes of cluster 
members (Icr-members as the most probable cluster stars, 2a- 
members as possible cluster stars, 3cr-members as possible field 
stars). Logically this classification is expressed by the following 
constraint: 1 A Pkin > Po A Pjk > Po A Pjh > Pq, 



where the threshold probability po =61' 
and 1 % for Sa-members. 



. for IfT-, 14% for 2a- 



P^P,■mm{Pk^n,PJH,PJK}, 



(6) 



3.4. Cluster parameter determination 

In general, the parameter determinatio n follows the rules applie d 
earlier to the clusters in the COCD ( Kharchenko et al.ll2005ah . 
However, due to the considerably fainter limiting magnitude, 
2MAst provides many more cluster members. In contrast to 
COCD-clusters we can, therefore, determine the cluster param- 
eters from the most probable members (Icr-members) located in 
the central cluster areas (r < ri, see Sect l3.4.ll for the defini- 
tion of ri), where the contamination by field stars is minimal. 
Consequently, the reliability of the results could be improved. 
Furthermore, the stars included in the solution had to fulfill the 
following conditions: 

- the 2MASS flag Qflg is "A" (i.e., signal-to-noise ratio 
S/N > 10) in each photometric band for stars fainter than 
Ks = 7.0; 

- the mean errors of proper motions are smaller than 10 mas/yr 
for stars with 5 > —30°, and smaller than 15 mas/yr for 
S < -30°. 

All parameters were determined in a process aiming at com- 
piling a consistent list of cluster members that in turn was used 
for revising the cluster parameters. This led to a series of itera- 
tions that each consisted of the following steps: 

1 . based on the parameters derived in the previous iteration, we 
computed for all stars in the cluster area the four membership 
probabilities discussed above (P,, Pkin, Pjh, Pjk)', 

2. based on the computed values of the probabilities, we con- 
structed a corrected list of cluster members and produced 
the spatial, kinematic, and photometric diagrams (shown in 
Figs. lCTTCll l: 

3. visual inspection and analysis of the diagrams for the de- 
termination of new cluster parameters: cluster coordinates 
and radii, average proper motion, age, distance, and redden- 
ing (as described in Sect. IIaTI 13131 and l3l3T l. If these 
parameters did not agree with those from the previous itera- 
tion and/or the member list differed from the previos one, the 
parameter list was corrected and the iteration was repeated, 
otherwise the process was stopped. 

As initial cluster parameters we used data taken from the litera- 
ture. If no literature data were available for a cluster, we set clus- 
ter parameters to values derived from an initial visual inspection 
of the data in the cluster field and to typical cluster parameters 
(ages, distances). 

To keep the derived parameters homogeneous and following 
the COCD practice, the pipeline was run by and under control of 
a single team member (NVK). A cluster was considered as "con- 
firmed" by the pipeline when compatible and non-contradictory 
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Figure 1. Comparison of apparent radii of open clusters (left 
panel for ri, right panel for derived with corresponding data 
from the COCD. The straight Une is a bisector, the dotted line is 
the Unear fit of the data. 



data on cluster structure, kinematics, and CMD were derived, 
and when the number of the most probable cluster members was 
found to be not less than five. Objects that did not fulfil these 
criteria were considered as "not confirmed". 



3.4.1. Cluster identification and spatial parameters of 
clusters 

Spatial cluster parameters were determined from stellar density 
profiles obtained from star counts. For each cluster, the counts 
were carried out in concentric circles around the cluster cen- 
tre with a minimum stepsize of 0?008. Fo r visual inspection 
and ch eck of cluster identification we used lArchinal & HvnesI 
(l2003j) as well as the SIMBAd^ and WEB DA databases, CDS 
PortaQ, and DeepSkyBrowseo resources. The density profiles 
were computed for four different stellar samples: Icr-, 2<t-, and 
3(T-cluster members, and all stars. The centres were determined 
by eye as the points of maximum surface density of the most 
probable cluster members (see Fig. lC.ll in Appendix ICl for illus- 
tration). Compared to the literature data, the coordinates of the 
cluster centres agree within 2.5 arcmin for more than 90% of 
the clusters. However, for the remaining clusters the differences 
can reach 0?1 and more. This is usually the case for relatively 
nearby large-sized clusters or associations, observed with a rich 
or irregular stellar background. For some objects this occurs due 
to mistakes in the source catalogues. 

The radial density profile (RDP) of the Icr-members was the 
decisive factor for the determination of the cluster size, though 
we also took into account their location in the cluster maps and 
the radial distributions of their /-Cs -magnitudes and proper mo- 
tions (see Fig. lC.ll in AppendixIClfor illustration). We assumed a 
centrally symmetric distribution of cluster members around the 
cluster centre and introduced three empirical structural param- 
eters To < ri < r2 describing the shape of the RDP, which 
were fitted by eye. The visible core radius tq corresponds to the 
distance from the cluster centre where the slope of the RDP be- 
comes flatter, the visible radius of the central part ri is the dis- 
tance where the decrease of the stellar density stops abruptly, 
and the actual (total) visible radius of a cluster r2 is defined as 
the distance from the cluster centre where the surface density of 
stars becomes equal to the average density of the surrounding 
field. To is a new feature considered here, while in our previ- 
ous work with the COCD we determined only ri and r2. Based 



Figure 2. Comparison of the average cluster proper motions 
(left panel for X-component, right panel for Y), derived in the 
present study with corresponding data from the literature. COCD 
clusters are shown by red circles, other clusters are shown by 
blue plusses. The straight line is a bisector, the dotted line is the 
linear fit of the data. 



on visual inspection of the corresponding data, tq, ri, and r2 
give empirical descriptions of cluster sizes, without a direct as- 
sociation to e.g. King's parameters Vc and rt (see Sect. 3.4.4.). 
However, we found that the relations tq < Vc < ri and rt > r2 
are generally obeyed. 

One can judge the accuracy of the derived radii by comparing 
our determinations with COCD-data. In Fig. [T] and Table |2] we 
compare the radii ri and r2 determined for MWSC and COCD 
clusters. One can see in Fig. [T]that both values are well repre- 
sented in logarithmic scale by a straight line 



p{COCD) =a + bp{MWCI) 



(7) 



http://simbad.u-strasbg.fr/simbad/ 
' http://cdsportal.u-strasbg.fr/ 
* http://messier45.com/ 



with fitting parameters as shown in Table |2] Although the 
MWSC radii are based on a much deeper catalogue (2MAst), 
they agree reasonably well with the results obtained with 
ASCC-2.5 data. 



3.4.2. Kinematic parameters 

For each cluster, the weighted mean components /J^ y of the 
cluster proper motion were computed from the proper mo- 
tions of the most probable members (Icr-members) within ri. 
Usually, we found 20-40 stars per cluster that fulfil this condi- 
tion. Compared to COCD, the typical formal accuracy of cluster 
proper motions is improved (0.7 mas/yr vs 0.4 mas/yr). No sys- 
tematic differences were found between literature and MWSC 
proper motions (see Fig.|2]and Table|2]i. 

We also compared MWSC proper m otions with dat a 
taken from the literature (as co llected in iDias et al ]j2002h . 
BasicaUy these are t he res ults of iDias et al.l (l200ll l2002h and 
Besheno v & Loktilil ( |2004|) . The cluster proper motions agree 
reasonably well. A higher dispersion can probably be explained 
by the fact that the data in the papers cited above are heteroge- 
neous compilations. 

As a rule, average radial velo cities of star clusters were 
adopted from Dias et al.l (l2002h . and lHarrisI (Il996l) . When possi- 
ble we computed mean radial velocities of the clusters using data 
on single stars from the CRVAD2-catalogue (iKharchenko et al.l 
[2007 ). Because there are few stars with known radial veloci- 
ties, we had to include data on Icr, 2cr and sometimes even 3cr- 
members distributed over the total area of a cluster within the 
radius r2- 
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Figures. Comparison of distance moduli (left 
panel), reddenings (middle panel), and ages 
(right panel) derived with literature data . 
Optical data dLoktin et alj|2001l : lLoktirJl2004l) 
are shown by red circles, the NIR data 
are marked by blue crosse s ([Fr oebrich et al 
K) and green plusses iBukowiecki et al 
20U|). The dotted lines labelled L(Loktin), 
F(Froebrich), and B(Bukowiecki) are the cor- 
responding fitted linear regressions (see Table[3] 
for their parameters). 



3.4.3. Photometric parameters and ages of clusters 

To determine the interstellar reddening, we used a colour-colour 
diagram {H — Kg), {J — H), while the diagram Qjhk^ i {J ~ 
Ks) was used for the control of the derived values. In the first 
diagram the reddening was estimated by a shift of the corre- 
sponding isochrone along the reddening vector from its intrin- 
sic position until the best agreement with the observations was 
achieved. In the second diagram we applied a horizontal shift 
of the isochrone which occurs when x is set to 2. Both diagrams 
are shown in Fig. lC.2l in Appendix|C]together with the initial and 
final positions of the isochrones. The distance to the cluster was 
determined by a shift of the isochrone in the colour-magnitude 
diagrams Kg, [J — H) and Ks^ {J — Kg). However, the choice 
of a proper isochrone presumes knowledge of the cluster age. 
Therefore, the cluster distance, age, and reddening were deter- 
mined in an iterative procedure. The iterations were stopped 
when the parameters did not change anymore. 

We found that the theoretical isochrones usually fitted the 
observed sequences in the Ks, (J — iiTs) -diagrams better than 
in Kg, (J — H). We observed a (J — H) colour shift that is 
typically of about a few hundreds of magnitude but it varied from 
cluster to cluster. We attribute this effect to a lack of a global 
calibration of the iJ-band in 2MASS (see Skrutskie et al...2006h 
and, additionally, to spatial variations of properties of absorbing 
dust that interfere in the iJ-band. Therefore, we introduce an 
empirical coiTection AH to provide a better isochrone fit in the 
Ks,{J — H) diagram. 

To determine the ages o f older clusters w e applied the proce- 
dure described in Kharche nko et al. I (l2005al) where the average 
age of the turn-off stars is used as a cluster age indicator. This ap- 
proach is justified for evolved clusters with well-defined turn-off 
points in the CMD (it applies for 33% of the considered clus- 
ters). In a younger cluster one observes a relatively poorly pop- 
ulated upper MS, which is augmented sometimes by a Pre-MS 
branch if the cluster is sufficiently close to the Sun. Therefore, 
the age determination of younger clusters was carried out with 
the isochrone-fitting technique. In both methods we used the 
Ks,{J ~ i^Ts) -diagrams built from the most probable members 



Table 2. Fitting parameters for relations between MWSC and 
literature data for cluster radii (in deg) and proper motions (in 
mas/y). 



Parameter 


N 


a 


b 


St.dev. 


logri 


126 


-0.23 ±0.04 


0.78 ± 0.04 


0.13 


log r2 


126 


0.01 ±0.03 


1.01 ± 0.04 


0.15 




189 


0.36 ±0.17 


0.87 ±0.04 


1.54 




189 


-0.23 ±0.19 


0.90 ±0.03 


1.92 



(Icr-members) located in the central cluster areas (r < ri). The 
selection of the best-fitting isochrone was made by eye. 

Comparing the newly derived distances and ages with COCD 
data, we find that both systems coincide for distances ( (Ad) — 
4 ± 25 pc) and show a small bias for ages ((Alogi) — 0.08 ± 
0.03). A comparison of the derived parameters with published 
data is sh own in Fig. [3] and Table [3j where we discuss the 
results of Loktin et al. (2001), Loktin (2004), Froebrich et ajJ 
(20 10), and Buko wiecki et al. (2011). The Loktin etal. {200^ 
an d lLoktinl (120041) results are based on photoelectric UBV obser- 
vations collected in WEBDA and provide the most extended and 
homogeneous set of cluster parameters obtained from optical 
photometry. Since our measurements and those of Loktin et al.l 
(2001) can be regarded as independent, we can consider the re- 
spective differences as a measure of accuracy of the derived val- 
ues. Assuming that the data have approximately equal weight, 
we find that our distances are accurate to within 16% (with a 
typical eiTor in the distance modulus of ±0.35 mag) and our 
reddening values E{B — V) are accurate to within ±0.06 mag 
(or about 10%). We note that for young clusters, which have 
almost vertical and poorly populated main sequences, the accu- 
racy of the par ameters is c onsiderably lower. The re sults from 
iFroebrich et all (120101) and iBukowiecki et al.l (l201lh were ob- 
tained with 2MASS data. Note that all these studies are based 
on a photometric selection of cluster members. 

To check possible systematic differences between our results 
and the published parameters, we constructed a linear regression 

Ap = p(lit) - p(MWSC) ^a + b [p(MWSC) - po] , (8) 

where p is one of the parameters (E{B — V), {V — My), or 
age log t), and po is a constant taken as 0.5 mag, 12 mag, and 8 
for these parameters respectively. Concerning the distance mod- 
uli and redd enings, the results generally agr ee reasonably well 
with optical (iLoktin et al.ll200ll:]Loktinl2004l) and with NIR data 

Table 3. Comparison of MWSC and literature data for photo- 
metric parameters and age. See text for the explanation of the 
parameters 



Parameter A'^ a b St.dev. Source 



iv~ 


Mv) 


109 


-0.06 ±0.05 


-0.03 ±0.05 


0.48 


L 


{V- 


Mv) 


43 


0.23 ±0.13 


-0.12 ± 0.13 


0.64 


F 


{V~ 


Mv) 


227 


0.42 ±0.08 


-0.38 ±0.08 


0.93 


B 


E{B 


~V) 


109 


0.00 ±0.01 


-0.01 ±0.01 


0.08 


L 


E{B 


-V) 


43 


0.02 ±0.05 


-0.17 ±0.05 


0.21 


F 


E{B 


-V) 


227 


-0.03 ±0.02 


-0.30 ±0.02 


0.16 


B 


log* 




109 


-0.15 ±0.03 


-0.18 ±0.03 


0.31 


L 


logi 




43 


0.55 ±0.09 


-0.54 ±0.09 


0.23 


F 


log* 




227 


-0.00 ±0.04 


-0.23 ± 0.04 


0.49 


B 



Refere nces: L=lLoktinl j2004h :F= lFroebrich et all 

i20T0l) ;B= lBukowiecki et alj ( 1201 ll) 
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log {r^ MWSC, deg) log (r, MWSC, deg) 

Figure 4. Comparison of tidal parameters of open clusters (left 
panel for Tc, right pa nel for rt), de rived with literature data. 
COCD cluster s ( Pisku nov et al ]|2007h are marked as red circles. 
iFroebrich et al.i 020071) data are shown by crosses. Their size and 
intensi ty is proportional to the quality flag in Froebrich et al. 
(l2007h ." Larger sizes and darker colours correspond to better 
quality. Two filled triangles show globular clusters identified in 
the 2nd quadrant. Their tidal parameters are taken for com pari- 
son from the 2010 edition of the catalogue o f lHartsl([T996l) . 

dFroebrich eta l.ll2010l) , but they show a systematic difference to 
the data of iBukowiec ki et al] (1201 lb . though the scatter of the 
data points is rather large. 

The comparison of ages points to a more serious problem 
that probably arises from the different approaches of member- 
ship determination. Indeed, a purely photometric selection usu- 
ally favours brighter stars as cluster members. However, they 
can be contaminating field stars, and their membership can be 
rejected by kinematic constraints. With such stars included as 
members, the cluster ages are underestimated, and the differ- 
ence to MWSC ages becomes negative, as seen in Fig. [3] for 
higher MWSC ages. On the other hand, when brighter mem- 
bers are excluded in the literature, the resulting age becomes too 
old compared to our estimate. Therefore, the membership ques- 
tion has a direct impact on the age determination of clusters, and 
this can be responsible for the discrepancies in the determination 
of {V — My) and E{B — V), too. Because our membership is 
based on kinematic and photometric constraints, we believe that 
our results are more reliable. Another reason for overestimated 
ages of young clusters can be the neglect of Pre-MS stars in the 
age determination (a procedure typical to some recent open clus- 
ter surveys), although they fix younger ages more reliably than 
MS-stars can do. As o ur analysis shows, all noticeable positive 
deviations of the B ukowiecki et al. clusters in Fig. [3] can 

be explained by the reasons above. 

3.4.4. Tidal parameters 

After completing the iterative parameter determination de- 
scribed in the previous subsections, we determined the tidal pa- 
rameters Tc (core radius), rt (tidal radius) using the same tech- 
nique as applied earlier to COCD clusters (see iPiskunov et al.l 
l2007l for details). The method is based on a three-parameter 
fit of cumulative King profiles to the observed density distribu- 
tion of cluster members. Since 2MAst is much deeper than the 
ASCC-2.5 catalogue, a considerably better member statistics in 
a MWSC cluster is now available. Already from this point of 
view we expect a better random accuracy of the tidal parame- 
ters. However, it is not a priori clear whether the new and old 
determinations do not differ systematically. 

In Fig. ID we compare the newly derived radii Tc and rt with 
our previous determinations (red open circles in Fig.|4|i. Though 



the data show a relatively large scatter for core radii Tc, we found 
no systematic differences between the new and old determinata- 
tions. Also, the tidal radii vt derived with ASCC-2.5 and 2MAst 
agree well for the majority of clusters. For a few nearby clusters 
(rt > 1 deg), the COCD tidal radii seem to be slightly overes- 
timated. However, we prefer to postpone the discussion on the 
significance of a possible bias until at least the next quadrant is 
processed. 

In Fig, m we also show the tidal parameters determined by 
IFroebrich et al.l(l2007b . who obtained cluster radii by fitting King 
profiles to the observed overdensities of 2MASS stars, neglect- 
ing any kinematic and/ or photom etric selection of members. The 
comparison with Froe brich et al is shown with crosses in 

Fig.|4] where the size and intensity of the crosses correspond to 
the quality level of the determinations as given by the authors. 
Larger crosses and darker colour describe a higher quality of the 
results. One can see that the clusters of Froeb rich et al. I (l2007l) 
exhibit systematically smaller angular core and tidal radii than 
the COCD clusters. This systematic displacement simply r eflects 
the fact that the clusters detected bv lFroebrich et al.l (|2007|) in the 
near-infrared are on average farther away from the Sun than the 
optical COCD clusters, whi ch represent the local cluster popula- 
tion. Considering only the Froebrich et all (120071) sample, we do 
not see a correlation of the results in both panels of Fig.|4] prob- 
ably due to the large relative errors. However, the higher disper- 
sion of the tidal radii of Froebrich et al. (2007) is in our opinion 
caused by a lo wer accuracy of their d etermination. This is not 
surprising since IFroebrich et aP (l2007[ ) fitted King curves to ob- 
served profiles contaminated by field stars. This introduced ad- 
ditional noise in their determination. Note that the clusters with 
the largest discrepancies in bot h core and tidal radii (F ig.lHl were 
marked as of lower quality bv lFroebrich et al.l ( l2007h . who also 
mentioned a generally much lower accuracy in their determina- 
tion of tidal radii compared to the core radii. 

In Fig. |4] we added comparisons for a few globular clus- 
ters shown as filled triangles. Out of four globular clusters lo- 
cated in the second quadrant we were able to fit density profiles 
for two objects (NGC 288 and Pal 2). The other two clusters 
(Pal 1 and Whiting 1) are too faint to show reasonable profiles 
in 2MAst. We compared our results with those fr om the cata- 
logue of the Galactic Globular cluster parameters (iHarrisll 19961 
2010 edition). The tidal parameters from the literature coincide 
sufficiently well with our determinations. 

To give the reader an overview of the properties of the open 
star clusters in the second quadrant, we show in Fig.|5]the distri- 
butions of the cluster parameters. 

4. Star cluster population in the 2nd quadrant of 
the Milky Way 

The current survey includes 652 objects (642 open, 2 globular 
clusters, and 8 associations) or 76% out of 871 from the input 
list, that could be confirmed as clusters in the sector of Galactic 
longitudes I ^ 90° . . . 180°. About half (122) of the rejected 
targets turned out to be random density enhancements of field 
stars, the others (97) are too faint to be detected in 2MAst. 

For each confirmed cluster we provide the basic set of clus- 
ter parameters (see Fig. |5] and Appendix). In Fig. |6] we show 
the distribution of the 652 clusters projected onto the Galactic 
plane {X, Y). For comparison, we also indicate the position of 
the grand design Perseus spiral arm with a pitch-angle of —6° fit- 
ted to the young COCD clusters by Piskunov et al. (2006). To es- 
timate the completeness of the MWSC cluster sample, we com- 
puted the surface density S of the clusters as a function of the 
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Figure 5. Distribution of star clusters in the second quadrant over derived fundamental parameters. 



projected distance from the Sun. The result is shown in Fig.|7]for 
650 open clusters and associations as well as for two subsamples 
of clusters, i.e. younger and older than logt = 7.9. Although 
the distribution does reveal fluctuations, it can be considered as 
flat up to dxY ~ 2.0 kpc. Then, the distribution drops, and the 
cluster density decreases slowly with lager distances. This can 
be interpreted as an indication of the completeness limit of our 
survey at dxY ~ 2.0 kpc. The completeness limit is defined 
by older clusters, while we observe an almost uniform distribu- 
tion of younger clusters in the range 2 kpc < dxY < 3 kpc, 
and a lower surface density only at larger distances. Probably, 
3 kpc is the completeness limit for young clusters in our survey. 
However, we cannot exclude that the decrease in the density be- 
yond about 3 kpc is a physical effect and related to the Perseus 
arm. A more conclusive answer is expected once the survey will 
be extended to ^ > 180°. 

Another interesting feature in the radial distribution of young 
clusters is a density peak at dxv ~ 0.6 kpc. This density in- 
crease can be seen also in Fig.|6]at I « 90° and is rather a foot- 
print of the star forming region in Cygnus than a mere density 
fluctuation. 



Considering the clusters within the completeness limit 
(dxY ~ 2 kpc), we found that the average surface density E 
is 94 kpc~^ in the 2nd quadrant, with contributions of 22 kpc~^ 
and 72 kpc^^ from clusters younger and older than log t = 7.9, 
respectively. This estimate is somewhat lower than the E = 
114 kpc~^, which we obtained w ith COCD clusters in a smaller 
area within dxY < 0.85 kpc (see lPiskunovetal.ll2006h . 



5. Summary 

We described the technique we applied to provide a global sur- 
vey of open clusters in the Milky Way and reported the first 
results achieved in the 2nd Galactic quadrant. The target list 
was compiled from present-day lists of known open clusters 
and cluster candidates but it also includes associations, mov- 
ing groups, and globular clusters. As the observational basis, 
we constructed a dedicated catalogue, called 2MAst, of about 
470 million stars with proper motions from PPMXL and pho- 
tometry from 2MASS. We developed a processing pipeline that 
provides simultaneous determination of the kinematic and pho- 
tometric membership as well as estimating cluster parameters. 
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Figure 6. Distribution of star clusters in the XF-plane centred at 
the position of the Sun. Blue squares mark the youngest (log t < 
7.9) clusters, black plusses show clusters with 7.9 < logi < 
8.3, black crosses correspond to clusters with 8.3 < logi < 
8.75, while red circles represent old clusters with \ogt > 8.75. 
The dashed magenta thick curve is the Perseus spiral arm fitted 
by a two-arm model of the spiral pattern with a pitch-angle of 
—6°. The large black dashed circle around the Sun indicates the 
completeness area of the COCD sample. The black solid circle 
marks the estimated probable completeness area of MWSC. 

The basic set of cluster parameters includes structure parame- 
ters (coordinates of the cluster centre, sizes), kinematic param- 
eters (average proper motions), photometric parameters (colour 
excesses, extinction in the Kg band, distance, and age). 

In total, 87 1 objects were examined in the 2nd quadrant (with 
850 cluster-like objects and 21 stellar associations). We con- 
firmed 652 of them (76%). Among them there are open and glob- 
ular clusters, cluster-like (compact) associations, moving groups 
and remnant clusters. We had to reject 12 large associations with 
sizes larger than a few degrees and with more than one peak 
in the density distribution of members over the projected sky 
area. For these "multi-centre" associations our pipeline failed to 
provide proper cluster parameters. Moreover, there were simply 
duplicate entries, or clusters too faint to be properly studied in 
2MAst. The majority of the rejected clusters, however, were ran- 
dom clusterings of field stars. 

The relatively high percentage of rejected clusters can be ex- 
plained both by the "risky" character of the input list, which 
includes all possible types of objects (including those that are 
treated as asterisms in the literature), and as a consequence of the 
enhanced control of the results, where the automated pipeline is 
extended by human decision at every stage. 

We determined the basic set of cluster parameters for all 652 
confirmed objects and were able to determine tidal parameters 
via a direct fit of King profiles to observed density distributions 
of cluster members. Radial velocities could be provided for 151 
objects (or 23% of the confirmed clusters). For the first time, 
distance and age estimates are given for 291 and 333 clusters, 
respectively. 
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Figure?. Distribution of the surface density S of star clusters 
in the 2nd Galactic quadrant versus their distance dxY from the 
Sun projected onto the Galactic plane. The density distribution 
of all clusters is given in black in both panels. In the upper panel 
we show the distribution of young clusters (logt < 7.9) in blue, 
in the lower panel that of older clusters(log t > 7.9) in red. The 
dotted lines indicate the completeness limit, the dashed horizon- 
tal lines correspond to the average density at dxY < 2 kpc for a 
given group of clusters. 
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Appendix A: MWSCS.Q2: The Catalogue of Milky 
Way Cluster Stars in the 2nd Galactic 
Quadrant 

The Catalogue of the Milky Way Cluster Stars in the 2nd 
Galactic Quadrant (MWSCS.Q2) exists in machine -readable 
form only and can be retrieved from the CDS online archiveQ 
The catalogue consists of 650 files with data on stars in areas 
with confirmed clusters. The information given for each star is 
shown in Table I A. fl 

Appendix B: MWSCD.Q2: the Catalogue of Milky 
Way Cluster Data in the 2nd Galactic Quadrant 

The Catalogue of the Milky Way Cluster Data in the 2nd 
Galactic Quadrant (MWSCD.Q2) exists in machine-readable 
form only and can be retrieved from the CDS online archive. 
The catalogue consists of the main table with the derived param- 
eters of the confirmed clusters, the list of all clusters, and a notes 
file. The latter two files contain all surveyed clusters including 
those that were NOT confirmed. To inform the reader on the data 
scope included in the catalogue, we describe here the main table 
(Table IrTT i. Object types are indicated by Al strings: empty - 
open cluster, a - association, g - globular cluster, m - moving 
group, n - nebulosity/presence of nebulosity, r - remnant clus- 
ter, * - asterism. Source object types are given with A3 strings 
describing a larger variety of previous classifications, e.g. dubi- 
ous clusters, clusters with variable extinction, embedded clus- 
ters, etc. The cluster metallicities have been adopted from the 

' ftp://cdsarc.u-strasbg.fr/pub/cats, [http://vizier.u-strasbgir| 



Table A.l. Contents of the MWSCS table of stellar data 



Col Label 


Units 


Explanations 


1 


RA 


h 


Right ascension J2000.0, epoch 2000.0 


2 


DE 


deg 


Declination J2000.0, epoch 2000.0 


3 


B 


mag 


B magnitude in Johnson system 


4 


V 


mag 


V magnitude in Johnson system 


5 


J 


mag 


J magnitude 2MASS 


6 


H 


mag 


H magnitude 2MASS 


7 


Ks 


mag 


Ks magnitude 2MASS 


8 


£j 


mag 


Error of J magnitude 


9 


£h 


mag 


Error of H magnitude 


10 


£Ks 


mag 


Error of Ks magnitude 


11 


PMx 


mas/yr Proper motion in RA*cos(DE) 


12 


PMy 


mas/yr Proper motion in DE 


13 


£PM 


mas/yr Error of proper motion 


14 


RV 


km/s 


Radial velocity 


15 


£RV 


km/s 


Error of RV 


16 


Qflg 




2MASS (ph_qual) JHKs photometric 








quality flag 


17 


Rflg 




2MASS (rd_flg) source of JHKs 








default mag 


18 


Bflg 




2MASS (bl_flg) JHKs components 








fit to source 


19 


2MASS 




2MASS (pts_key) unique source 








identifier in catalogue 2MASS 


20 


ASCC 




ASCC-2.5 number 


21 


Sp 




Spectral type and luminosity class 


22 


Rci 


deg 


Angular distance from the adopted 








cluster centre 


23 


Ps 




Spatial probability of cluster 








membership 


24 




% 


Proper motion probability of cluster 








membership 


25 




% 


Photometric JKs probability of cluster 








membership 


26 


PjH 


% 


Photometric JH probability of cluster 








membership 


27 


No 




MWSC number 



compilation of lDias etal ] (l2Q02h . their uncertainties are of about 
0.1 dex. 

Appendix C: MWSCA.Q2: the Atlas of Milky Way 
Cluster Diagrams in the 2nd Galactic 
Quadrant 

For each cluster we prepared two pages of diagrams that are 
combined in the atlas. For a better understanding, we describe 
the diagrams of the atlas taking the cluster King 19 as an exam- 
ple. 

We consider the cluster map as the basic diagram of the first 
page, and the CMDs Ks,{J — H) and Kg, (J — Kg) as the basic 
diagrams of the second page. Stars are shown as coloured circles 
or dots. Symbols and their colours have the same meaning in 
all plots. Cyan symbols mark stars outside the cluster radius r2, 
green symbols stars within r2 - The most probable kinematic and 
photometric members (Icr-members) are indicated in black for 
members located within ri, red for members between ri and r2, 
and blue for stars outside r2- Cyan bars show the uncertainty for 
Icr-members (page 2). 

Page 1 of the atlas (Fig. IC.ll ) contains five diagrams with 
spatial information and a legend on the derived cluster parame- 
ters. The right panel is a map of the cluster surrounding, the left 
panels show magnitudes Kg, proper motions PM^, PMy, and 
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Table B.l. Contents of the MWSCD main table 



rnl T ahel 




Units 


XJiAUId.lId.llUl is 


1 No 




— 


MWSC number 


2 Name 




— 


NGC, IC or other common designation 


3 type 






object type 


4 RA 




h 


RA J2000.0 of the centre 


5 DE 




deg 


Dec J2000.0 of the centre 


6 1 




deg 


Galactic longitude of the centre 


7 b 




deg 


Galactic latitude of the centre 


8 ro 




deg 


Angular radius of the core 


9 ri 




deg 


Angular radius of the central part 


10 r2 




deg 


Angular radius of the cluster 


11 PMx 




mas/yr 


average proper motion in 
RA*cos(DE) 


12 PMy 




mas/yr 


average proper motion in DE 


13 epM 




mas/yr Error of PM 


14 RV 




Ifln/s 


Average radial velocity 


15 erv 




km/s 


Error of RV 


16 nRV 




— 


Number of stars used for RV 


17 Nlarg 






Number of most probable (Icr) members 
in ro 


18 Nlari 




— 


Number of most probable (Icr) members 
in ri 


19 A^lcTra 






Number of most probable (Icr) members 
in r2 


20 d 




pc 


Distance from the Sun 


21 E{B~ 




mag 


Colour-excess in [B — V) 


22 (Tfs - 


Mk,, 


) mag 


Apparent distance modulus 


23 E{J - 


Ks) 


mag 


Colour-excess in (J — Ks) 


24 £;( J - 


H) 


mag 


Colour-excess in (J — H) 


25 AH 




mag 


AH 


26 logf 




logyr 


Logarithm of average age 


27 eiogt 




logyr 


Error of log t 


28 Nt 






Number of stars used 
for the calculation of log t 


Zy Tc 




pc 


King' core radius 


30 




pc 


Error of core radius 


31 rt 




pc 


King' tidal radius 


32 £,.t 




pc 


Error of tidal radius 


33 fc 




pc-^ 


King' normalisation factor 


34 Ek 




pc-2 


Error of normalisation factor 


35 source 






Source for MWSC list and input 
parameters 


36 source type 




Source object type 


37 [Fe/H] 






Metallicity 



responding radius; weighted average components PMx,y of 
proper motion with their rras errors and number of stars used to 
compute the average; the average radial velocity, RV, rms error, 
and the number of stars used to compute the average; distance to 
the cluster, d, distance modulus, {Kg — Mk,)', NIR interstellar 
reddening, E{J — H), E{J — Ks), and interstellar extinction, 
A{Ks)', cluster age, its rms error, the number in brackets gives 
the number of stars used to compute the average age, or it is - 1 
if an isochrone fitting was applied. AH shown below the pho- 
tometric diagrams indicates the empirical correction to the H- 
magnitude introduced in Sect. l3.2l and l3.4.3l 

Page 2 (Fig. lC.2l l contains three diagrams with kinematic in- 
formation (left panels), and six diagrams with photometric infor- 
mation (right panels). 

The three left panels with kinematic data: the two upper dia- 
grams show PMx,Y vs. Ks relations, i.e.'TM-magnitude equa- 
tion". Magenta vertical lines correspond to the average proper 
motion of the cluster. The magenta dashed line shows the ap- 
parent magnitude if™'^, which coiTesponds to the bluest colour 
( J — i^s) of the adopted isochrone. The bottom panel is the vec- 
tor point diagram of proper motions. 

The six right panels with photometric data: the two upper di- 
agrams are CMDs {Kg, [J—H] and Kg, {J —Kg))- The magenta 
curve is the apparent isochrone closest to the determined clus- 
ter age. Solid blue lines outline a domain of 100% photometric 
members. Solid red lines are the ZAMS and TAMS (shown only 
in Ks,{J ~ Kg)). The magenta dashed line shows the apparent 
magnitude K™''. The thic k yellow circles mark the s tars used for 
the age determination (see lKharchenko et al.ll2005al for details). 
The black aiTows show the vectors of increasing extinction. The 
four bottom panels show the two-colour {H — Kg)/ [J — H) di- 
agram (left column) and Qj^f^-colour diagram (right column). 
The upper row is for stars brighter than K™'^, the lower row is 
for stars fainter than K"''^-. Magenta curves indicate the apparent 
isochrone (i.e., observed colours), while cyan curves show the 
intrinsic isochrone. 

The legend is the same as in page 1 . 



surface density N versus distances r of stars from the cluster 
centre. 

The sky map: stars are shown by circles. The size corre- 
sponds to star brightness arranged in six Kg magnitude bins. The 
blue cross indicates the cluster centre determined in this study, 
while the blue plus sign is the cluster position taken from the 
literature. If by chance other clusters appear in this cluster area, 
their centres are marked by magenta plus signs. Large blue cir- 
cles (shown by dotted, solid, or dashed curves) indicate the clus- 
ter radii rp, ri, or r2, respectively. The left panels: blue vertical 
lines (dotted, solid, or dashed) mark ro, ri, or r2. Magenta hori- 
zontal lines in the PM vs. r diagrams correspond to the derived 
average proper motion of the cluster The RDPs in the bottom 
panel are shown with green for all stars, blue for Su-members, 
magenta for 2cr-members and black for Icr-members. 

The legend gives cluster name, MWSC number, and COCD 
number in parentheses; equatorial RAj2ovio^ Decj2aoa, and 
galactic /, b coordinates of the cluster centre; apparent cluster 
sizes ro, ri, r2 and number of Icr-members within the cor- 
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Figure C.l. Open cluster King 19 in the MWSC atlas (page 1). See text for the description. 
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Figure C.2. Open cluster King 19 in the MWSC atlas (page 2). See text for the description 
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